ASYMPTOTICS OF COEFFICIENTS OF MULTIVARIATE 
GENERATING FUNCTIONS: 
IMPROVEMENTS FOR SMOOTH POINTS 
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Abstract. Let X],3eN d Fp x ^ be a multivariate power series. For example ^ Fpx 13 could 
be a generating function for a combinatorial class. Assume that in a neighbourhood 
of the origin this scries represents a noncntirc function F = G/H p where G and H 
are holomorphic and p is a positive integer. Given a direction a 6 N 1 ^ for which the 
asymptotics are controlled by a smooth point of the singular variety H = 0, we compute 
the asymptotics of F na as n — > oo. We do this via multivariate singularity analysis 
and give an explicit formula for the full asymptotic expansion. This improves on earlier 
work of R. Pemantle and the second author and allows for more accurate numerical 
approximation, as demonstrated by our examples. 



I. Introduction 

Let 'Yjj3m d Fp %13 be a multivariate power series. For example Fpx 13 could be a gener- 
ating function for a combinatorial class. In |PW02t IPW04] Pemantle and Wilson derived 
asymptotic expansions for the coefficients Fp as /3 — > oo for large classes of series that 
arise often in applications. In this article we further their program of asymptotics of 
coefficients of multivariate generating functions. 

Assume that in a neighbourhood of the origin the power series Ylf3m d Fp x>3 is the 
Maclaurin series of a nonentire function F = G/H p where G and H are holomorphic and 
p is a positive integer. For example F could be a rational function. Using multivariate 
singularity analysis we derive the asymptotics of F na for a G and n — > oo in the case 
that these asymptotics are controlled by smooth points of the singular variety H = 0. 

Our presentation is organized as follows. In Section [2] we set our notation and basic 
definitions. In Section [3] we give two explicit formulas for all the terms in the asymptotic 
expansion of F na . This improves upon |PW02] which gave formulas for the leading term 
only and only for the case p — 1. Furthermore, we prove that the expansions for F na 
are uniform in a as claimed in [PW02J. In Section H] we go on to express our formulas in 
terms of the original data G and H for easier use in calculation. In Section we apply 
our formulas to three combinatorial examples and demonstrate the greater numerical 
accuracy of many-term asymptotic expansions over single-term expansions. Lastly, in 
Section [6] we present the two theorems from analysis on Fourier-Laplace integrals that we 
use in our proofs. 

2. Notation and definitions 

Let N + and 1R + denote the set of positive natural numbers and positive real numbers, 
respectively. For r G M and k G N = N + U {0} set r k = r(r + 1) • • • (r + k — 1), the kth 
rising factorial power of r, with the convention that r° = 1. For m G N + , x G C m , and 
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i < m let Xi denote component i of x and x = (xi, . . . ,x m _i). For a G N+, x G C m , 
and n G N + define a + 1 = (ai + 1, . . . , a m + 1), a! = • • • a m \, na = (nai, . . . , na m ), 
x a = x" 1 ■ ■ ■ x^ , and d 01 = d" 1 ■ ■ ■ <9,^ m , where dj is partial differentiation with respect 
to component j. For c G C m let -D(c) = {x G C m : Vj < m \xj\ < \cj\} and C(c) = 
{s G C' m : Vj < m = |cj|}, the polydisc and polycircle centred at the origin with 
polyradius (|ci|, . . . , |c m |), respectively. Finally, all sum indices will begin at the zero 
tuple unless indicated otherwise. 

Fix d G N + and let G, H : T> — ► C be holomorphic functions on a polydisc PCCf 
Assume that G and H are relatively prime in the ring of holomorphic functions on T>. Let 
p G N + and define F = G/H p . Then F is holomorphic on T> \ V, where V is the analytic 
variety {x G P : = 0}. By |Sch05t Example 4. 1.5, Corollary 4.2.2] the variety V has 

dimension d — 1 and X> \ V is a domain, that is, a nonempty open connected set. Assume 
G £> \ V and V ^ 0. Let J^pem F P X& be the Maclaurin series of F, so d^F(0)//3\ = Fp 
for all /3 G N d . 

We will derive asymptotics for Fp as (5 — > oo along straight lines through the origin 
and off the axes, that is, for F na with a G and n — * oo. For d = 2 and /3 — >• oo along 
more general paths see }Lla06] . 

First we recall some key definitions from [P W02j . Just as in the univariate case, the 
asymptotics for the Maclaurin coefficients of F are determined by the location and type 
of singularities of F, that is, by the geometry of V. Generally the singularities closest to 
the origin are the most important. We define 'closest to the origin' in terms of polydiscs. 

Definition 2.1. Let c G V. We say that c is minimal if there is no point x G V such 
that \xj\ < \cj\ for all j < d. We say that c is strictly minimal if there is a unique 
x G V such that \xj\ < \cj\ for all j, namely x = c, and we say that c is finitely minimal 
if there are finitely many such values of x. 

The variety V always contains minimal points. To see this let c G V and define 
/ : V fl -D(c) — > R by f(x) = \fx\ + ■ ■ ■ + x 2 d . Since / is a continuous function on a 
compact space, it has a minimum, and that minimum is a minimal point of V. 

We focus on the singularities of F with the simplest geometry, namely the regular or 
smooth points of V. For a summary of what is known about non-smooth points see the 
survey |PWj . 

Definition 2.2. A point c G V is called smooth if V is a smooth complex manifold in 
a neighbourhood of c. 

Equivalently, a point c G V is smooth iff djH(c) ^ for some j iff Cj is a zero of order 
one of z i — H(ci, . . . , Cj_i, z, Cj + i, . . . , q) for some j. See [BK861 page 363] for a proof of 
the first equivalence. 

We will approximate F na with integrals and, in so doing, will need to consider the 
singularities of F relevant to a. These singularities are called critical points. For a more 
geometric explanation of their relevance and the stratified Morse theory behind it see 
[PW] . 

Definition 2.3. Let a G and c G V be a smooth point. We say that c is critical 

for a if it is a solution of the system of d — 1 equations 



a 1 1 x\diH(x) = . . . = a d l XdddH{x). 
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We say c is isolated for a if, in addition, it has a neighbourhood in which it is the only 
critical point for a. 



When H is a polynomial, the system of d equations in d unknowns given by H(x) = 
and the critical point equations generally has a finite set of solutions. 

Remark 2.4. At times it will be convenient to work in projective space. Recall that 
CP d_1 is the set of equivalences classes of C d \ {0} under the equivalence relation ~ given 
by x ~ x' iff x = Xx' for some A G C \ {0}. Let ~~ : \ {0} — > CP d_1 be the natural map 
which takes a point x to its equivalence class [x]. 

While the definition of F na only makes sense for a G N d , the definition of critical point 
is well-defined for all a G CP d ~ l with at least one nonzero component. 

3. The full asymptotic expansion 

Pemantle and Wilson [PW02] showed that if a G Nl and c G V is strictly minimal, 
smooth, critical and isolated for a, and nondegenerate (which we will define shortly), 
then there exist G C such that for all N G N + one has the asymptotic expansion 



P — r -na 
- 1 na 



b k n-^' 2 ~ k + O [ n -^' 2 - N ) 

lk<N 

as n — > oo. They also derived a similar expansion for degenerate points in the case d = 2 
and gave an explicit formula for 6 f° r a U d when p = 1. 

In this section we derive explicit formulas for all bk and all p. To formulate our results 
we employ the following functions. 

Definition 3.1. Let c G V be smooth, and assume without loss of generality that 
ddH(c) 7^ 0. By the implicit function theorem there exists a bounded neighborhood W 
of c and a holomorphic function h on W such that (w, h{w)) G V, ddH(w, h{w)) ^ for 
all w G W, and h(c) = c^. Suppose q ^ 0, so that we may also assume that h is nonzero 
on W. _ 

For j < p define Uj : W -»• C, E : [-1, l]**" 1 -> C^ 1 , and u,-, h : S" 1 ^ R C(c)) -> C 

by 

ttj = Uj o E 
h = hoE. 

Furthermore, for c critical for a define g : — > C by 

5 W = iog(|4) +i 2=*.. 

Then ^ is well-defined for any a G CP d_1 with 7^ since 5^ does not depend on the 
magnitude of a and since h is nonzero on W. Moreover, Uj, h, and g are all C°° functions. 
If det <f' (0) 7^ 0, then c is called nondegenerate. 

In the context of a single pair of appropriate c and a we will use the functions of 
Definition 13. II without further introduction. 
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Uj (w) = lim (-y)-P+J— ((y-h(w)yF(w,y) 



Now for the first theorem. 



Theorem 3.2. Let d > 2 and a G N+. If c G V is strictly minimal, smooth with 
c d d d H(c) 7^ 0, critical and isolated for a, and nondegenerate, then for all N G N, 

-1/2 



(*) 



i<p-i fc<iv 



(p- 1 - 



+ O ( n P-l-(d-l)/2-") 



as n — > oo. 
Here 



H l+k {u$){t ) 



K2k 



l) k 2 l+k l\(l + k)V 



g(t) = g(t)-g(t )-±(t-t )g"(t)(t-t ) T , U is the differential operator -J2 r ,s(9"(t ) l ) r , s d r d s , 
and to — 0. In every term of L k (uj, 'g) the total number of derivatives of w and of g" is at 
most 2k. 

In the case d = 2 we can drop the nondegeneracy hypothesis. 

Theorem 3.3. Let d = 2 and a G Nl. If c G V is strictly minimal, smooth with 
CdddH(c) 7^ 0, critical and isolated for a, and v > 2 is least such that <7^(0) ^ 0, then 
for all iVeN, 



(*0) 



- 1 na 



(aa^n) 



EE 



(a d n + 1) 



p-i-j 



7TC ^ ' ^ ' (/' ' 1 - /)! /! 

j"<P fc<7V V ^ J 7 J 



■(a rf n)-^ir(«^) 



O ( n p-i-(2AM-i)A>) 



as n — > oo for t> even and 

'(Ho^n)- 1 ^ 



(41) 



27rv ^— ' , z -f r (p — 1 — j)!j! 

j<p fc<iV v ^ Jy J 



+ 0(n p - 1 -( iV+1 ^) 



as n — > oo for v odd. 
Here 



2fc+t>Z+l' 



K2fe 



:-i) f r(^ 

/!(2fc + vZ)! 



dt 



2k+vl 



{ujg l )(t ), 



_ / y\l-p( k+vl+i \ 



Kfc 



Z!(ife + uZ)! 



I |a| 1/,; isgna— ) (wj£')(*o), 



fc+uZ 



g{t) = g{t) - g(t ) - a(t - t ) v , a = g^(t )/v\, ( = e in/( - 2v \ and t = 0. In every term of 
L^ en (u,g) the total number of derivatives of u and g^ is at most 2k, and in L^, dd (u,g) 
at most k. 

In the formulas above we take z~ 1 ' v = \z\~ 1 l v e~ ia,rgz / v for z £ C \ {0} with argz e 
[-7r/2,7r/2]. 

To prove Theorems 13.21 and 13.31 we follow the same general approach as in [PW02J and 
summarized in the following steps: (1) use Cauchy's integral formula and strict minimality 
to express c na F na as a <i-variate contour integral over a contour almost touching c; (2) 
expand the contour across c and use Cauchy's residue theorem along with the smoothness 
of c to express the innermost integral as a residue; (3) calculate the residue explicitly, and 
take the resulting (d — l)-variate contour integral and change to real coordinates to get a 
Fourier-Laplace integral; (4) use theorems from analysis (see Section [6]) to approximate 
the integral asymptotically. 

Lemma 3.4 ( \PW02\ proof of Lemma 4.1], Steps 1 & 2). Let a e N d + and c e V. If 

c is strictly minimal and smooth with CdddH(c) ^ 0, then there exists e G (0, 1) and a 
polydisc neighborhood D of c such that 

F na = (2^J x ^ld W + 0(en 

as n — > oo, where X = D fl C(c) and is the residue of y h- > F(w,y)y~ adn ~ 1 at 

Lemma 3.5 (Step 3). In the previous lemma, 



Thus 



c™F na = (2 7 r) 1 - d ]T (p^tjlli / uj(t)e- a ^dt + 0(e"), 
j<p ^ 3)-3- Jx 

as n — > oo, where X = i? _1 (X). 

Proof. This is a straightforward calculation. Let w E X . Since c is smooth, h(w) is a 
simple zero of y t— > i?(u», y) and so a pole of order p of ?/ t— > -F(w, ?/). Since /i is nonzero 
on iy, /i(w) is also a pole of order p of y \— > F(w, y)y~° ldn ~ 1 . Thus 

fl « = t^w (IT' ( fa - ^w*- ^'^i 



v ~ i\ A ^ V 



from which the first identity follows by definition of Uj. Thus 

Jx ™ na+1 



j<?-i (p _ 1 _ ^ w "° v ft(c) ' n 



Km«2 Wr ' 



(via the change of variables w = E(t)) 



which with Lemma [3.41 proves the stated formula for c na F na . □ 

Remark 3.6. In the case d = 1 Lemma 13.51 simplifies: h, Uj, and R become 0-ary 
functions, that is, constants (h becomes c), and there is no integral. Thus we arrive at 
the following known result. 

If c G V is strictly minimal and smooth, then there exists e G (0, 1) such that 



as n — > oo, where 



lim(-a;)- p+J (J^-J ((x - c) p F(x)). 



Moreover, if c G V is finitely minimal and smooth and every point of V D C (c) is smooth, 
then the asymptotic expansion of F n is the sum of the expansions around each point of 

vnc(c). 

Before proceeding to Step 4 we collect a few technical facts about g. 

Lemma 3.7. Let a G N+ and c G V. If c is strictly minimal, smooth with CdddH(c) ^ 0, 
and critical and isolated for a, then for all t G X we have dfg(t) > with equality only 
at t = and g'(t) = iff t = 0. 

Proof. Firstly, g(0) = by definition. Secondly, 9?g(t) > with equality only at 
since c is strictly minimal. Lastly, by the implicit function theorem, d m h(w) = 
—d m H(w, h(w))/ddH(w, h(w)) for all m < d and w G W. So for t G X we have 

B a(t)--\r r it m 1 d m H(E(t),h(E(t))) a m 
m9{) ~ m h(E(t)) d d H(E(t), h(E(t))) + V 

Therefore g'{t) = iff 

a" 1 c m e iim d m H{E{t), h(E(t))) = h(E(t)) d d H(E(t), h(E(t))) 

for all m < d iff {E(t), h(E(t))) is critical for a iff t = since c is isolated for a. □ 



Proof of Theorem \3.S\ (Step 4)- By Lemmas 13.41 and 13.51 there exists e G (0, 1) and an 
open bounded neighbourhood X of such that 



= (2tt) ^ T^T—SuT 1 ^ + ( e 



as n — > oo, where Jj in = j^Uj{t)e~ adn ^dt. 

Choose k G C^°(X) such that k = 1 on a neighbourhood Y of 0. Then 

/ ifB = / K (t)u j {t)e- a * ,iS ®dt + I {I- K{t))uj{t)e- adn ^dt. 
J X Jx 

The second integral decreases exponentially as n — > oo since 9ftg is strictly positive on 

the compact set X \ Y by Lemma 13.71 By Lemma 13.71 again and our nondegeneracy 
hypothesis, we we may apply Theorem 16.11 with to = to the first integral. Noting 
that Lk(htUj,g) = Lk(uj,g) because the derivatives are evaluated at and k = 1 in a 
neighbourhood of 0, this gives 

^ ^ ' ' k<N 

= {^y id ~ 1)/2 (detg"(0)y 1/2 J2(a d n)- k L k (u v g) + 0(n~^- N ) 

k<N 

as n — > oo. Hence 

n™F foi-JV (a d n + iy- l -\ , „« 

c F na =(27r) 2^ (p-i- J ^ + °( e ) 
j<p yF J, ' J ' 



E E ^T~T[T (J 2 ™^" 1 det ^(°)) " 1/2 (««t")"%(5i, 9) 



(p — 1 — 

j<p k<N v ^ J ' J 

p-l-(d-l)/2-N^ 



+ o( 



n 



as n — > oo, as desired. □ 



The proof of Theorem 13.31 is similar but uses Theorem 16.21 instead of Theorem 16.11 
In the case of finitely minimal points of V we simply take an open set W for each 
finitely minimal point so that W contains no other finitely minimal points, repeat the 
proofs above for each such W, and add the asymptotic expansions. Thus we have the 
following. 

Corollary 3.8. Let a G and c G V. If c is finitely minimal and every point 
of V PI C(c) satisfies the hypotheses (excluding strict minimality) of Theorem 13.21 or 
Theorem 13.3} then the asymptotic expansion of F na equals the sum of the expansions 
around each point of VnC(c). □ 

Finally, we show that the asymptotic formulas for F na are uniform in a. This was 
claimed without proof in |PW02] . 

Proposition 3.9. Let d > 2 and K C MP d_1 be compact. Suppose that for all a G K 
there exists a unique c G V that is strictly minimal, smooth with c d d d H(c) ^ 0, critical 
and isolated for a, and nondegenerate. Suppose that all these points c lie in a bounded 
open set V C V such that x d d d H(x) ^ for all x G V. Then for each N G N the big-oh 
constant of (Jt) stays bounded as a G N+ varies and a stays in K. 
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Proposition 3.10. Let d = 2 and if C RP 1 be compact. Suppose there exists v > 2 
such that for all a G X there exists a unique c G V that is strictly minimal, smooth with 
CdddH(c) 7^ 0, critical and isolated for a, and t> is the least integer greater than 2 such 
that g^(0) 7^ 0. Suppose that all these points c lie in a bounded open set V C V such 
that XdddH(x) ^ for all x G V. Then for each JVeN the big-oh constant of (4*0) and 
(4*1) stays bounded as a G N + varies and a stays in if. 

Proof of Proposition ^. We take up where the proof of Theorem 13.21 left off. Recall that 
g depends on a. Let us emphasize this dependence by writing g a . By Theorem 16.11 it 
suffices to show that for any fixed positive integer P there exists M > such that for all 
a G if and all (3 G N d_1 with \(3\ < P we have 

(<?) ||^ Q ||oo = sup{|^ Q (t)| : t G X} < M, 

where X C R d_1 is a suitable neighborhood of 0. 

To prove such a bound we first show that the correspondence between a direction in 
if and its critical point is continuous. To this end it will be helpful to introduce the 
logarithmic Gauss map 7 : V _> CF d - x defined by x 1— > [xidiH(x), . . . ,XdddH(x)}. Note 
that 7 is well-defined since (xxdiH(x), . . . , XdddH(x)) 7^ on V by hypothesis. Note also 
that a G if is critical for c iff 7(c) = a. 

Since ii is holomorphic on V, each o^-ii is holomorphic on V, so that 7 is continuous. 
By hypothesis, for every a G if the preimage 7 _1 ({a}) contains exactly one element. 
Thus 7~ : : if — > V is well-defined. Moreover 7 _1 is continuous since 7 restricted to 
7~ 1 (if) is a continuous bijection on a compact space into a Hausdorff space and therefore 
a homeomorphism. 

Each a G if has associated to it a point 7~ 1 (a) G 7 _1 (if) which has has associated 
to it the functions of Definition 13.11 We now show that finitely many such functions 
will do to handle all points of 7 _1 (if). For each c G 7 _1 (if) let D be the polydisc from 
Lemma 13.41 containing c and h the nonzero holomorphic function on D associated to c. 
Let D' be a polydisc such that D' C U' C D. The collection of all such D' x h(D') 
forms an open cover for the compact space 7 _1 (if). So this cover has a finite subcover 
D[ x hi(D[), . . . , D[ x hi(D[). Let Di, . . . , A be the superpolydiscs corresponding to 
D[, ...,£>;. To handle the various E of Definition EH define e : Ylj<iDj x h 1 , l] d_1 
(£<z-i ( W) £) i_> (u' 1 e 1 * 1 , . . . , Wrf.xe 1 * 1 *- 1 ). Then e is continuous, and for all j the open set 
e~ l (Dj) contains D'- x {0}. Since D'- is compact, by the tube lemma [Mun75l Lemma 
5.81 there exists a neig hborhood Yj of in [-1, l]^ 1 such that e-\Dj) D x Yj. Set 

Y = f]j<i Y j C [-1, l] d_1 , and let X be a neighborhood of with X C X CY. Thus to 

each a G if are associated some c = 7~ 1 (a), _D', /ij, ft,^, and log ( ) +i ^ m<d — £ m — 

g» a , where the last two functions are defined on Yj and so on Y. 

With this setup we now show ( ( v ) ). Since Aj : Y — > M for j < / and 5 : if x Y — > R 

defined by A^i) = logfrT^r) an d B(a,t) = ^2 m<d —t m are continuous and X 
and if are compact, we have that Ma '■= Y2j<i max{A, (t) : t G X} < 00, Mb '■— 
max{i?(o;,t) : a G if, t G X} < 00, and 

I |<7a| loo < sup Aj(t) + B(a, t) (for some j < I) < M A + M B , 

a£K,t£X 



la, oo 



a bound that is independent of a. Similarly, since each hj is C°° over Y, each \\d^g 
for < P stays bounded as a varies within K, as desired. □ 
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The proof of Theorem 13. 101 is similar but uses Theorem 16.21 instead of Theorem 16.11 



4. Rewriting the expansion in terms of the original data 

To actually compute with the formulas in Theorems 13.21 and 13.31 it is helpful to rewrite 
the quantities involved in terms of the original data G and H . The propositions below 
give formulas for calculating g"(0) and Uj in terms of derivatives of G and H. 

Proposition 4.1 ([RWa, Theorem 3.3]). Let a G and c G V. If c is smooth with 
CdddH(c) 7^ and critical for a, then for all l,m < d with I^mwe have 

9 (0)lm 



c 2 d (d d Hf 



x (dmHdiH + caiddHdmdiH - d m Hd d diH - d^d^H + d 2 d H\ 

9"(0)u =^|| + ^jj- 2 (m? + c d (d d Hd?H - 2d lH d d d lH + { -^fdjH)) 
c d o d H c d (d d H) z V o d H J 

□ 

In the presence of symmetry Proposition 14.11 simplifies greatly. 

Proposition 4.2 ([RVVal Proposition 3.4]). Let a G and c G V. If x i— > H(x) is 
symmetric, a has all of its components equal, and c lies in the positive orthant, is strictly 
minimal, smooth with c d d d H(c) ^ 0, and critical for a, then c has all of its components 
equal and for alH, m < d with I ^ m, 



g"(0)lm = q, g"(0)u = 2q, and det <?'(()) = dq 



d-l 



where q = l + ^(djH - d^H) 



□ 



Proposition 4.3. If c G V is smooth with c d d d H(c), then for all j < p and to G W we 

have 

Wj(iy) = {—h(w)) p J hm 



y-*h(w) dyiQ(w,y)P 1 

H(w,y) 
y-h(w) ' 

lim dP d Q{w,y) = -^-<^ +i #(w, h(w)). 

y-+h{w) K -+- I 



where Q : ((W x C) fl V) \ V — > C is given by Q(w,y) = ^j™ ■ Moreover, 

1 .-,+■! 



In particular, 



«o(wj = 7 ^j. 

— h(w)d d H(w, h(w)) 



Proof. The first statement is just the definition of Uj. 
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Let (w, y) e ((W X C) n V) \ V. Then H(w, h(w)) = and 




H(w,y) y d%H(w,h(w)) 



(y-h(w)) 



n 




Therefore Q extends to a holomorphic function on (W x C) fl T> and 




(y-h(w)) 



□ 



5. Examples 



Let us apply the results of Section [3] and H] to several combinatorial examples. We used 
Maple 11 to do the calculations, and our worksheets are available at http : //www . cs . auckland . ac . nz/~rai 
First we mention two shortcuts to finding strictly minimal points. 

Proposition 5.1 ([PW., Theorem 3.16]). If the coefficients of F are all nonnegative 
and there is a critical point for a, then there is a minimal critical point for a in Mi.. □ 

A rf-variate power series a a x a is called aperiodic if the Z-span of {a G N d : a a ^ 0} 
equals Z d . 

Proposition 5.2 ([PW, Proposition 3.17]). If 1 — H is aperiodic and has nonnegative 
coefficients, then every minimal point of V is strictly minimal and lies in R^. □ 

Example 5.3 (d = 2, p = 1, N = 2). Consider the bivariate generating function 



F(x 1 ,x 2 ) = - 



1 — X\ — X 2 — X\%2 



whose coefficients -P/3 1)( g 2 are called Delannoy numbers and count the number of lattice 
paths from (0, 0) to (fa, j3 2 ) with allowable steps (1, 0), (0, 1) and (1, 1). We compute the 
first two terms of the asymptotic expansion of F na as n — > oo for a = (3, 2). 



The critical points of V are 



(-| + |Vl3,-| + iVl3) and (-§ - |^13, -§ - |^13). 



Both points are smooth and the first point, which we denote by c, is strictly minimal by 



Propositions 15.21 and 15.11 



Applying the results of Sections [3] and HJ we get 



F na = (cr 3 c 2 - 2 ) n (b n-^ + hn- 3 /' + O (n^)) 

as n -f oo, where q 3 c 2 2 w 71.16220050, b = + VlS) « 0.3690602772 and 

h = - ^Ifj^l (79^ + 767) w -0.01853610557. 
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Comparing this approximation with the actual values of F na for small n (using 10-digit 
floating-point arithmetic), we get the following table. 



n 


1 


2 


4 


8 


16 




25 


1289 


4.673345- 10 b 


8.527550909-10 i:i 


3.978000114-10^ 


c- na (fe n" 1/2 ) 


26.26314145 


1321.542224 


4.732218447-10 6 


8.581184952-10 13 


3.990499094-10 28 




24.94407138 


1288.354900 


4.672799360-10 6 


8.527311037-10 13 


3.977972633-10 28 


one-term relative error 


-0.05052565800 


-0.02524610085 


-0.01259771042 


-0.006289501355 


-0.003142026054 


two-term relative error 


0.002237144800 


0.0005004654771 


0.0001167557713 


0.00002812906104 


0.000006908245151 



For an arbitrary a, the two-term asymptotic expansion of F na is just as easy to compute 
symbolically in a. The corresponding constants Ci, C2, Oi, a<i are square roots of rational 
functions of aci, a^, and \Ja\ + of. The exact formulas are somewhat long, so we omit 
them. 



Example 5.4 (d = 

generating function 



2, p = 1, N = 5). Fix q G (0,1) and consider the bivariate 

F(x ll x 2 ) = - ■ — 

1 — qx\ + qx\%2 — xfx2 

which arises in the context of quantum random walks. Motivated by an example from 

[BP], we compute the asymptotics of F na as n — > oo for a = (2, 1 — q). 

There is one critical point of V for a, namely, c := (1,1). This point is smooth. 

Explicitly solving for x 2 as a function of X\ in H(xi,X2) = and applying the minimum 

modulus theorem shows that c is minimal. However, c is not finitely minimal, because 

for every pi G C(l) there exists P2 G C(l) such that (pi,P2) G V. By a modification of 

Lemma [3.51 described in [BP} Theorem 3.2] we may still apply Theorem 13.31 to c. Doing 

so, we get 



oo, where b 



as n 

n^ 1 , and n _4//3 terms are zero. 



= &o^" V3 4 

(l-9) 2/3 
32/3 g l/3(l +g )l/3r(2/3) 



6 4 n- 5/3 + 0(n- 2 
and 6 4 = 



) 



3 1 / 6 (g 4 +22g 2 +l)r(2/3) 
280g 5 /3( g+ i)5/3( 1 „ 9 )2/3 7r - 



The n- 2 / 3 , 



Comparing this approximation with the actual values of F na for small n (using 10-digit 
floating-point arithmetic), we get the following table for q = 1/2. 



n 


2 


4 


8 


16 


32 


bon- 1 / 3 +b 4 n~ 5 / 3 
one-term relative error 
two-term relative error 


0.1875000000 

0.1953794677 

0.1855814246 

-0.04202382773 

0.01023240213 


0.1523437500 

0.1550727862 

0.1519865960 

-0.01791367352 

0.002344395487 


0.1221771240 

0.1230813520 

0.1221092630 

-0.007400959937 

0.0005554313097 


0.09739671811 

0.09768973380 

0.09738354495 

-0.003008476011 

0.0001352526066 


0.07744253816 

0.07753639314 

0.07743994970 

-0.001211930578 

0.00003342426606 



Example 5.5 (d = 3, p < 3, N = 2). Consider the (d + l)-variate generating function 

W(xi, ...,x d ,y) = ^-r, 

l-yB(x) 

where A(x) = l/(l - J2j=i^+x), B i x ) = 1 - (1 - ex(x))A(x), and ei(x) = J2i=j x j- 
Using the symbolic method (as presented in [FS, Chapter 1], say) it is not difficult to 
show that W counts all words over a <i-ary alphabet A, where Xj marks occurrences of 
letter j of A and y marks occurrences of snaps, non-overlapping pairs of duplicate letters 
counted from left to right. Here A(x) counts snapless words over A, the so-called Smirnov 
words. For more details see [RWbj. 
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The coefficient W( n ,...,n,s) is then the number of words with n occurrences of each letter 
and s snaps. For n G N + let ip n be the random variable taking a word over A with n 
occurrences of each letter and returning the number of snaps in the word. We compute 
the expectation and variance of ip n as n —>■ oo. If a — (1, ... , 1), then 



(l-ei(x)H 
/ ) 



d 2 W 
dy 2 



(x,l) + ^(x,l)) 



A(x)- 2 5(x)(5(x) + 1)(1 - ei(x))- 3 



'(1 - e 1 (x))~A 

V / 7 

E(va-E(^ n ) 2 . 



Let if (x) = 1 — ei(x). Then the only critical point of V is c := (1/d, . . . ,1/d), and 
it is strictly minimal by Propositions 15.11 and 15.21 Applying the results of Sections [3] 
and S to := W(x,l) (with p = 1), F 2 {x) := dW/dy(x,l) (with p = 2), and 

F 3 (x) := <9 2 Wy<9?/ 2 (x, 1) + dW/dy(x, 1) (with p = 3) with d — 3, we obtain 



E(W, 



E(^ 



[(2 ? rn) 2 rf]- 1 / 2 E J<2 E fc<2 ^g^^L fc (n 2 ,^) + g 
[(2nnyd)]-^E 3<1 E k< 2 



ii 



-2\ 



(-i)i' 



h L k (u 2 ,j,g) + 0(n- 3 ) 



yj n -i _ Y^- 2 + o( n -3) 

i^-l + oK 1 ), _ 

[(2 7 rn) 2 d]-V 2Ej<3Efc<2 (l^ 



k L k (u 3J ,g) + 



n 



n 



-3> 



\{2^nfd}- V2 E K1 £ fc<2 ^ggf n-*L fc (5x + O 

|Jn- 2 + 0(n-3) 
^n 2 -fn + 0(l), and 
|n + 0(l). 



Comparing these approximations with the actual values for small n (using 10-digit 
floating-point arithmetic), we get the following table. 
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n 


2 4 8 


(3/4)n 

(3/4)n - 15/32 
one-term relative error 
two-term relative error 


1.000000000 2.509090909 5.520560294 
1.500000000 3 6 
1.031250000 2.531250000 5.531250000 
0.5000000000 0.1956521740 0.08684620409 
0.03125000000 0.008831521776 0.001936344398 


(9/16)n 2 

(9/16)n 2 - (27/64)n 
one-term relative error 
two-term relative error 


1.800000000 7.496103896 32.79620569 
2.250000000 9 36 
1.406250000 7.312500000 32.62500000 
0.2500000000 0.2006237006 0.09768795635 
0.2187500000 0.02449324323 0.005220289555 


(9/32)n 
relative error 


0.800000000 1.200566706 2.31961973 
0.5625000000 1.125000000 2.250000000 
0.2968750000 0.06294253008 0.03001342380 



6. Asymptotics of Fourier-Laplace integrals 

We conclude with the main analytic results behind the proofs of Theorems 13.21 and 
13.31 Pemantle and Wilson's approach in [PW02] and our approach here to deriving 
asymptotics for the Maclaurin coefficients of F requires an asymptotic expansion for 
integrals of the form 

/ u(t)e- W9{t) dt 
Jx 

as uj —>■ oo, where X C R d ~ l is open and u and g are complex valued functions on X. 
Pemantle and Wilson used Watson's lemma and Morse's lemma to prove the existence 
of a full asymptotic expansion for these Fourier-Laplace integrals as they are commonly 
called but gave an explicit formula for the leading term only. In contrast we use the two 
theorems below which give explicit formulas for all terms. 
All function spaces mentioned are complex valued. 

The first theorem deals with stationary and nondegenerate points, that is, points to 
such that g'(to) = and det(g"(to)) 7^ 0, respectively. 

Theorem 6.1 ( |Hor83l Theorem 7.7.5]). Let X C R^ 1 be open, u £ C 6 °°(X), and 

g £ C°°(X). If 9ftg > 0, $lg(t ) = 0, g has a unique stationary point t £ suppu, and t 
in nondegenerate, then for every N £ N + there exist M > such that 

I u(t)e-^dt = e-"*™ (det \ ^ ^ uj~ k L k { U) g) + O (uj-^I 2 -*) 

Jx V V 71 J J k<N 

for uj > 0. Here L k is the function defined in Theorem 13.21 (but without the stipulation 
to = 0). Moreover, the big-oh constant is bounded when the partial derivatives of g up 
to order 3(iV + \(d — 1)/2~|) + 1 all stay bounded in supremum norm over X. □ 

The second theorem deals with degenerate stationary points for d = 2 and is proved 
by adapting Hormander's approach. 

Theorem 6.2 ((10 Theorem 1]). Let X C R be open, u £ C C °°(X), and g £ C°°(X). 
If Ifog > 0, $tg{to) = 0, g has a unique stationary point to £ suppn, and v > 2 is least 
such that g( v '(to) ^ 0, then for every N £ N + there exists M > such that 

( u(t)e-^dt = 2 ( auj y /V uJ- 2k/v Lr n {u, g) + O {u^ 2N+1 ^ v ) 

^ X V k<N 
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for io > and v even, and 



/ u(t)e-^dt = e-^ito) ^) 1/V u~ k ' v L odd {u, g) + (u 

^ X V k<N 



-(N+l)/v 



) 



for <jj > and v odd. Here L^ en and L^ dd are the functions defined in Theorem 13.31 (but 
without the stipulation t = 0). Moreover, the big-oh constants are bounded when the 
derivatives of g up to order (v + l)(iV + 1) + 1 all stay bounded in supremum norm over 
X. □ 
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